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Abstract. Random magnetic flares above the accretion disks of Active Galactic Nuclei can account for the pro- 
duction of the primary radiation and for the rapid X-ray variability that have been frequently observed in these 
objects. The primary component is partly reprocessed in the disk atmosphere, forming a hot spot underneath the 
flare source and giving rise to distinct spectral features. Extending the work of Czerny et al. (2004), we model the 
fractional variability amplitude due to distributions of hot spots co-orbiting on the accretion disk around a super- 
massive black hole. We compare our results to the observed fractional variability spectrum of the Seyfert galaxy 
MCG-6-30-15. According to defined radial distributions, our code samples random positions for the hot spots 
across the disk. The local spot emission is computed as reprocessed radiation coming from a compact primary 
source above the disk. The structure of the hot spot and the anisotropy of the re-emission are taken into account. 
We compute the fractional variability spectra expected from such spot ensembles and investigate dependencies 
on the parameters describing the radial spot distribution. We consider the fractional variability Fvar with respect 
to the spectral mean and the so-called point-to-point definition Fpp. Our method includes relativistic corrections 
due to the curved space-time in the vicinity of a rotating supermassive black hole at the disk center; the black 
hole's angular momentum is a free parameter and is subject to the fitting procedure. We confirm that the rms- 
variability spectra involve intrinsic randomness at a significant level when the number of flares appearing during 
the total observation time is too small. Furthermore, the fractional variability expressed by -FVar is not always 
compatible with Fpp. In the special case of MCG-6-30-15, we can reproduce the short-timescale variability and 
model the suppressed variability in the energy range of the Ka line without any need to postulate reprocessing 
farther away from the center. The presence of the dip in the variability spectrum requires an increasing rate of 
energy production by the flares toward the center of the disk. The depth of the feature is well represented only 
if we assume a fast rotation of the central black hole and allow for considerable suppression of the primary flare 
emission. The modeled line remains consistent with the measured equivalent width of the iron Ka line complex. 
The model can reproduce the frequently observed suppression of the variability in the spectral range around 6.5 
keV, thereby setting constraints on the black hole spin and on the disk inclination. 
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1. Introduction 

The variability of the X-ray emission is one of the basic 
characteristics of Active Galactic Nuclei (AGN) . Well es- 
tablished since the EXOSAT observations (Pounds 1985), 
the variability is expected to be an important tool in 



^ Email: agata@camk.edu.pl 

" Email: anne-marie.dumont@obspm.fr 



Goosmann et al.: Structure and X-ray spectra of AGN accretion disks III: rms- variability 



the study of the accretion process, although its aperi- 
odic noise-like character was early on recognized as a 
source of interpretation difficulty (McHardy & Czerny 
1987; Lawrence at al. 1987). 

The phenomenological approach to X-ray variations 
has brought significant information on the properties of 
the lightcurves, such as their power spectral density (e.g. 
Markowitz et al. 2003; McHardy et al. 2005) and their 
process nonlinearity (Uttley et al. 2005). Nevertheless, the 
physical processes behind the observed variability have not 
yet been fully understood. The same is true for accreting 
galactic black holes. The main source of the ambiguity 
is that several different models can account for the av- 
erage radiation spectra, e.g. the 'hot inner flow model' 
(Ichimaru 1977; Narayan & Yi 1994; Kato et al. 2004; 
Sobolewska et al. 2004a), the 'continuous corona model' 
(Liang & Price 1977; Bisnovatyi-Kogan & Blinnikov 1977; 
Czerny & Elvis 1987; Haardt & Maraschi 1991; Czerny et 
al. 2003; Merloni 2003), the 'patchy corona model' (Galeev 
et al. 1979; Haardt et al. 1994; Stern et al. 1995; Ross et al. 
1999; Nayakshin et al. 2000; Sobolewska et al. 2004b), or 
the 'jet-basis model' (Henri & Pelletier 1991; Markoff et al. 
2005). Several models can also account for the power den- 
sity spectra, e.g. the 'shot model' (Lehto 1989), the 'self- 
organized critical state model' (Mineshige et al. 1994), 
or the 'traveling perturbations model' (Lyubarskii 1997; 
King et al. 2004) . Apparently, it is impossible to uniquely 
determine the accretion flow geometry on the basis of 
power spectra or from the analysis of the mean spectral 
shape alone. 

Therefore, a 2-D approach is needed, which would com- 
bine the time and the energy dependence. Such an ap- 
proach was pioneered in the context of binary black holes 
using various techniques, e.g., phase lags (Miyamoto et 
al. 1988; Cui et al. 1997), correlation functions (Nolan 
et al. 1981; Poutanen & Fabian 1999), coherence func- 
tions (Vaughan & Nowak 1997), frequency resolved spec- 
tra (Gilfanov et al. 2000; Zycki 2002; Zycki 2003), wavelet 
analysis (Lachowicz & Czerny 2005). Now the appropriate 
data is starting to be available for AGN (e.g. Yaqoob et 
al. 2003; Kaastra et al. 2004; Ponti et al. 2004; Markowitz 
2005). We therefore concentrate on modeling the energy- 
dependent fractional variability amplitude (Vaughan et 
al. 2003). Long light curves obtained from Chandra and 
XMM-Newton satellites allow us to determine this func- 
tion for a few AGN. Such detailed information provides a 
challenge to the models. 

In the present paper, we consider the magnetic flare 
model of the X-ray emission, suggested by Galeev et al. 
(1979) in the context of Cyg X-1, and further developed in 
various directions by several authors, both in the context 
of galactic sources (Poutanen & Fabian 1999; Beloborodov 
1999) and AGN (Haardt et al. 1994; di Matteo 1998; 
Merloni & Fabian 2001; Torricelli-Ciamponi et al. 2005). 
In this picture, the magneto-rotational instability (MRI) 
operating in the disk occasionally leads to the rise of large 
magnetic loops high above the disk surface. Field recon- 
nection occurs in this region and leads to the creation 



of localized regions of hot plasma (flares). The flares are 
cooled by Comptonizing the disk emission, which pro- 
duces a power law X-ray radiation, customarily called pri- 
mary emission. A fluctuating number of flares depending 
on various parameters is supposed to be responsible for 
the observed AGN variability. Recent 3-D simulations of 
the magneto-rotational instability within the disk seem to 
support such a picture (e.g. Miller & Stone 2000; Turner 
2004) , but simulations including the Compton cooling are 
still to be done. 

The hard X-ray emission of a flare leads to the cre- 
ation of a hot spot at the disk surface where the X-ray 
emission is reflected/reprocessed. The combined emission 
of all spots thus constitutes the so-called reflected compo- 
nent in the total X-ray spectrum, and the issue of the re- 
processing of the incident X-ray radiation by the disk sur- 
face has been extensively studied (e.g. Lightman & White 
1988; Ross & Fabian 1993; Zycki et al. 1994; Nayakshin 
et al. 2000; Ballantyne et al. 2001; Rozahska et al. 2002). 

Here, we extend our previous work on this subject 
(Czerny et al. 2004). We improve our description of the 
spot emission by taking its dependence on the distance 
from the spot center and on the local emission direction 
into account. We use these new computations to model the 
standard fractional variability amplitude as well as the 
point-to-point fractional variability amplitude. We com- 
pare the results to observational data of MCG-6-30-15. 

2. The model 

2.1. Radiative transfer within tine spot 

Let the flare be located at a certain height above an ac- 
cretion disk. It illuminates the disk surface. We assume 
that the flare source is very compact and approximate it 
by an isotropic, point-like source. The emitted spectrum 
is of a power law shape, with the photon index F = 1.9, 
extending from 1 eV to 100 keV. 

The reprocessed emission depends on the distance from 
the flare. Therefore, we divide the patch into concentric 
rings (see Fig.0. Within each ring, the incident angle of 
the hard X-ray illumination and irradiating flux is approx- 
imately constant. It is convenient not to parameterize the 
rings by their radial distance from the projected position 
of the flare (spot-center), but instead by /i, which is the 
cosine of the incident angle 6 of the irradiating flux (see 
Fig. nj. We vary jjl — cos 9 in constant steps of 0.1, ensur- 
ing that all rings intercept the same fraction of the total 
illuminating flux. 

The reflected spectrum of each ring is calculated for 
five values of the local emission angle, ip, measured with 
respect to the disk axis. We choose equal steps in cos'0 
from 0.1 to 0.9. The local radiation spectrum is computed 
accurately by performing radiative transfer simulations in 
plane-parallel approximation. The vertical structure of the 
illuminated disk is computed using the code of Rozaiiska et 
al. (2002). Following an argument made by Collin et al. 
(2003), we assume that the total flare duration is signif- 
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Source of 
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Fig. 1. Spot approximated by five rings below a flare. In 
each of the rings, the reflected spectrum is calculated for 
the corresponding incident angle 9 of the X-ray radiation 
and the ionization parameter at the surface. 

icantly shorter than the characteristic timescale for the 
restoration of the hydrostatic equilibrium in the suddenly 
irradiated disk. We thus use the vertical density profile of 
the disk as computed before the flare turns on. 

The local radiative transfer computations are con- 
ducted by coupling the codes titan and noar, as de- 
scribed in Dumont ct al. (2000) and updated in Dumont 
et al. (2003). To model the accretion disk atmosphere, 
we assume a vertically stratified, plane-parallel medium 
being irradiated by hard X-ray radiation from above its 
upper surface. From underneath, low-temperature black- 
body radiation is added to account for the heating of 
the underlying disk. The computations are performed for 
a 10®Mq black hole with a disk accretion rate of 0.001 
of the Eddington accretion rate. A distance of 18 Rg 
{Rg = GM/c^) from the black hole is assumed, and the 
ratio of the incident flux to the disk flux equals 144 at 
the center of the spot. Details of these calculations are 
described in Goosmann et al. (2006). 

The overall response of the irradiated material to the 
incident flux is in agreement with known results: a hot 
layer forms at the top of the disk atmosphere, that is 
roughly at the inverse Compton temperature, followed by 
a steep transition to colder, less ionized layers. Such be- 
havior is due to the requirement of the hydrostatic equi- 
hbrium (Field 1965; Krolik et al. 1981; Raymond 1993; 
Ko & Kallman 1994; Rozahska & Czerny 1996; Nayakshin 
et al. 2000; Ballantyne et al. 2001; Rozahska et al. 2002). 
However, the horizontal stratification of the spot intro- 
duced in our model brings new aspects to the resulting 
reprocessing: we include the fact that the disk surface is 



most strongly irradiated directly below the flare, while the 
irradiating flux decreases and the incident angle of the ir- 
radiation increases with the distance from the hot spot 
center. Therefore, the surface inside a hot spot shows a 
significant radial gradient of its properties, including the 
locally emitted spectrum. The incident angle and the local 
incident flux are different in each ring. Large differences 
between the ring spectra are seen in the soft X-ray band. 
Spectra from the outer rings are much harder than those 
from the inner rings, since the reflector there is less ion- 
ized. The high-energy parts of the local spectra from all 
rings are rather similar to each other. 



2.2. Ring integration and dependence on the local 
emission angle 

The final spectrum of a spot, for a given local emission 
direction, is determined by the integration over the rings. 
Due to the adopted parameterization of the rings, the in- 
tegration actually reduces to a simple summation. If the 
fiare emission is somehow collimated toward the disk, the 
outer parts of the hot spot are not irradiated. Therefore, 
we can treat the opening angle of the irradiating cone as 
an additional parameter of the model. Here, we adopt a 
half-opening angle 9o = 60°. This value of Qq corresponds 
to half of the illuminating flux being reprocessed by the 
disk medium and also ensures that our calculations re- 
main accurate with respect to the local conditions of the 
assumed hydrostatic equilibrium. The resulting local spec- 
tra are shown in Fig. El These V'-dependent spot spectra 
are further adopted as representative of computations of 
a whole accretion disk covered by many spots. 

We note that by integrating over the spot surface, we 
average the spectrum and lose some spatial resolution. 
We do that in the present model for simplicity. In real- 
ity, the properties of the incident and reprocessed radia- 
tion change across the spot, which affects the single spot 
lightcurve. An investigation of this effect is part of a de- 
tailed study about individual flare spots (Goosmann et al. 
2006). From this study, it follows that in our case the in- 
tegration over the spot surface does not induce a major 
error. The intensity of the local spectra at a given emis- 
sion angle ip changes regularly across the spot surface and 
only within a factor of ^^ 2 — 3. Also, the spectral slope 
only changes by a factor < 2. 

The shape of the radiation spectrum of a whole spot 
does not depend strongly on the emission angle ip, while 
the overall normalization does. The normalization rises 
toward higher values of ip, which represents a limb- 
brightening effect. The spectrum for all emission angles are 
much harder than the local spot spectrum used by Czerny 
et al. (2004). This is caused by our present assumption of 
short lasting flares, while Czerny et al. (2004) assumed 
that the flare duration is long enough for the underlying 
irradiated disk layers to settle into a new hydrostatic equi- 
librium. 



ms- variability 




Fig. 2. The spot spectra (intensity) integrated over the 
rings for 5 values of the emission angle ip: cosTp — 0.1, 0.3, 
0.5, 0.7, and 0.9. The upper spectrum is the side view. 
Parameters adopted: M — lO^M©, R ~ 18 Rg, riidisk — 
0.001, and Lx/id = 144. 



2.3. Flare distribution across the disk 

To define the flare distribution across the disk, we gen- 
erally follow the simple approach outlined in Czerny et 
al. (2004). We assume that the spectral shape of all flares 
is the same, since it consumes too much computer time 
to repeat such complex computations separately for spots 
located at various disk radii. On the other hand, we take 
the fact that the reflected (spot) spectrum shows strong 
anisotropy of the emissivity (limb-brightening effect) into 
account. 

The source is characterized by rimcam the mean num- 
ber of all coexisting flares at a given time. The spots are 
assumed to appear randomly across the disk surface, be- 
tween the inner radius Rm (equal to the marginally stable 
orbit around a black hole) and the adopted outer radius, 
i?out- The probability distribution of their radial and az- 
imuthal position, {Ri,(j)i), is determined by the general 
law 



p(i?i,0i)=p(i?i)p(0i) = 



(7,.ad + 2)R^- 



1+1 






+2 



RZ 



^2 "" 2^- 



(1) 



This distribution reduces to the uniform distribution 
over the disk surface, if 7rad = 0, while positive (negative) 
values of 7i.ad lead to an enhanced number of flares in the 
outer (inner) disk parts. 

The spot radius, i?x, is assumed to be the same for all 
spots. The flare luminosity and, consequently, the local 



spot radiation flux generally scales with the flare position, 
i?i, as 



^inr — ^n 



Rir 



-Aa 



(2) 



The dissipation within the Keplerian disk scales with 
R^^, so the typical value of /3rad is expected to be of that 
order. 

The duration of the flare may also depend on the flare 
location, so we describe it as 



^lifc — ^lifcn 



-Rin 



(3) 



Here, (5rad = means that the flare lifetime is indepen- 
dent of the flare location, while for S^^d — 1.5, the flare 
lifetime scales with the local Keplerian timescale. 

The observational appearance of the source also de- 
pends on the duration of the observation (i.e., integration 
time) and on the inclination angle of an observer with re- 
spect to the disk symmetry axis. The value of the black 
hole mass, M, is another relevant parameter since all dis- 
tances and timescales scale with M. 

The flare distribution is fully described by the model 

parameters: Af, Tlmoan, -Rin, ^out, Rx, Fq, /3rad, 7rad, <5rad, 

and tiifoQ. The spot radius i?x can be replaced, if conve- 
nient, by the mean X-ray luminosity of the source, Lx- 

2.4. Modeling sequences of spectra 

The observed properties of the spot distribution are 
strongly affected by special and general relativistic effects. 
These effects must be taken into account with the maxi- 
mum possible accuracy since they account for a significant 
fraction of the resulting variability. Even a stationary spot 
on a Keplerian orbit gives a complex variability curve if 
the relativistic effects are included. Light bending causes 
a magnification of the emission from a relatively small 
disk area (located roughly behind the black hole). Flares 
passing through this area strongly contribute to the rapid 
fiuctuations of the signal. 

We apply the code KY (see Dovciak et al. 2004 and 
Dovciak 2004 for a description) here to compute the rel- 
ativistic corrections between the local spot spectra and 
the actually 'observed' ones. The code incorporates all 
aspects of the light propagation from between the disk 
surface and a distant observer. It includes the effects of 
light bending, as well as gravitational and Doppler shift- 
ing of the photon energy. The ruling parameters for the 
relativistic effects are the viewing angle i of a distant ob- 
server, measured with respect to the disk normal, and the 
dimensionless Kerr parameter a of the black hole. The 
anisotropic character of the local (TA-dependent) spot spec- 
tra (see Sect. I2.2|) is also relevant and included in our 
modeling. 

The resulting disk spectrum seen at a given viewing an- 
gle i is based on a combination of various angle-dependent 
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local spot spectra. However, the integration over the disk 
surface for a large number of hot spots is computation- 
ally time-consuming, so that we cannot produce detailed 
lightcurves for all spots. However, to model the data, such 
lightcurves are actually not needed because the spectra 
can only be measured in bins Tobs of the typical order of 
a few thousands of seconds. Therefore, within each time 
bin, we consider the appearance/disappearance of the hot 
spots and their motion. Then we construct the emissiv- 
ity 'belts', which represent the traces of the spots on the 
disk surface over their duration. This method is identical 
to the one applied in Czerny et al. (2004). Finally, the 
spectrum 5*1; (-E) of the radiation received by a distant ob- 
server within a time bin k, is computed. We calculate N 
of such sequences for a given flare distribution. Usually 
we adopt TV = 50 since this corresponds well to the best 
observed lightcurves covering a few hundreds of thousands 
of seconds. 

Another modification, introduced in the current model 
in comparison to the method of Czerny et al. (2004), is the 
survival of flares between consecutive time bins. Czerny et 
al. (2004) treated each time bin as being separate from the 
previous one, and all spot positions were newly generated. 
This is correct if the time bins do not come from a single, 
uninterrupted observation or if the effective life time of 
the flares is shorter than the length of a time bin. In the 
current modeling, we follow the positions of the hot spots 
over consecutive time bins. A randomly sampled position, 
according to an adopted spot distribution, is assigned only 
to new born flares. This Monte-Carlo sampling also con- 
siders an adopted birth rate, which is determined by the 
mean number of flares, their life time, and by the spatial 
spot distribution. This parameterization allows an appro- 
priate description of the overall X-ray properties of the 
source. 

2.5. Fractional variability 

The obtained radiation spectra sequences allow us to cal- 
culate the fractional variability for a given flare distribu- 
tion. We follow the definition used by Ponti et al. (2004): 



N 



.(^) = E 



k=l 






(4) 



where N is the number of time bins, S'k(-B) is the inte- 
grated radiation spectrum for a single time bin fc, and 
{S{E)) is the mean spectrum. We also determine the sta- 
tistical errors for this determination of i^var , as it is done 
for the observed data (see Ponti et al. 2004). Of course 
the modeled lightcurve is not limited by the S/N ratio, 
although it is limited by the adopted total duration of the 
observation and by the width of the energy channels for 
the final -FVar. 

By statistical errors, we mean the error determination 
from a single simulated lightcurve. In the simulations, we 
use a fine energy grid of 800 points in a 2 - 10 keV energy 
range. The dimensionless Fvai(£^) is determined at each 



energy point. When we compare the model to the data, 
we bin the results using broader energy bins. If n original 
energy points contribute to the new energy bin centered 
around -Ek, we introduce 



-Fvar(-Ek) 



Z]"=l-P'var(-Bi) 



(5) 



and the error of this quantity is determined as 



(5Fvar(-Bk) = 



T:UFUEi)~nFUE^) 



iFl,{E^) + l)FUE^) 
n{N - 1) 



1/2 



(6) 



The first term in this expression determines the error of 
the mean value by averaging over the energy bins and as- 
suming an error of Fvar (Ei) — Fvar (-^k ) for a given i^var (Ei). 
The second term is the expected statistical error for a 
normalized variance effectively calculated from n{N — 1) 
points. Both components are necessary, since in the model 
the values in nearby energy bins are not statistically in- 
dependent. 

2.6. Point-to-point fractional variability 

We also calculate the point-to-point fractional variability 
(Edelson et al. 2002; Vaughan & Fabian 2004), which em- 
phasizes the character of the shortest timescale variation. 
This quantity is calculated by 



p2 ,j^. ^ y-' 4[^k+i(i^)-^k(i^)] 
h (iV-l)[5k+i(i?) + 5k( 



S^.{E)Y 



together with the rebinned version 



Fvp{E\d 



e: 



-.iEpp(.Ei) 



(7) 



(8) 



and the appropriate statistical errors SFpp. They are 
computed from the same formula as (Si^var, replacing 
Fvar(S) by Fpp{E) in Eqn. Eland i^var(^) by Fpp{E) re- 
spectively. Our definition of Fpp is even more local than 
the one previously used because it adopts the local mean 
instead of a global mean at each step. 



3. Results 

In our model setup, the vertical disk structure and the lo- 
cal spectra, as a function of the emission angle ip, are ac- 
curately computed for a distance of 18 Rg from the black 
hole. These local spectra are then applied to hot spots 
across the whole disk. On the one hand, this is an impor- 
tant simplification in our modeling. On the other hand, 
however, it helps to reduce the complexity of the problem. 
The model setup described in the previous section focuses 
on investigating the effects on the variability caused by 
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General Relativity, by the radial spot distribution across 
the disk, and by the anisotropic local spectrum emitted 
by the spots. The shape of the local spectra for each spot 
are thus assumed to be independent of the spot location, 
while their normalization is ruled by the radial emission 
profile of the disk. Further radial dependencies, like the 
power of the primary emission illuminating the disk or 
the hydrostatic equilibrium of the disk atmosphere, are to 
be considered in future work. 

The following presentation of our results starts with 
a comparison to the previous approach of Czerny et al. 
(2004) in Sect. O Then, in Sect. ES we point out the 
main effects on the variability as a function of the model 
parameters, and we explain them qualitatively. Finally, in 
Sect. 13.31 we apply the model to some variability data of 
the Seyfert-1 galaxy MCG-6-30-15. 

3.1. The effect of the spot structure and of the 
anisotropic spot emission 

Since we made two important modifications of the model 
in comparison to the approach of Czerny et al. (2004), we 
show the significance of these changes in Fig. |3| for two 
sets of parameters - one corresponding to a Schwarzschild 
black hole and one corresponding to a maximally rotating 
Kerr black hole. The mean number n„iean of flares is ad- 
justed so that for both choices of the black hole spin, the 
variability level is roughly the same (see Sect. I3.2|l . We 
compare the results obtained from (a) the code of Czerny 
et al. (2004), (b) the new code taking into account the 
cloud survival between the time bins, but with the old spot 
radiation spectrum, and (c) the new code, considering the 
spot structure and the angle-dependent spot emission. 

The differences between cases (a) and (b) are mainly 
in the normalization. The survival of the flares between 
the consecutive time bins enhances the variability since 
it decreases the randomness in the cloud position. The 
discrepancy between the two approaches increases with a 
rising ratio of the flare life time iufe to the bin size Tobs- 
The considered examples were calculated for tufc being 
independent from the disk radius (^rad =0). If the radial 
dependence is strong (e.g. 5rad — 1-5), the cloud survival 
also induces some difference in the shape of Fvar . 

The most profound difference is introduced by the 
change of the local spectrum. Case (c) shows a much 
stronger energy dependence because the new local reflec- 
tion spectrum has more prominent spectral features, in- 
cluding the iron Ka line (see Sect. 12.21) . This is promising 
in the perspective of reproducing the data. The model of 
Czerny et al. (2004) predicted a far too shallow energy 
dependence, in comparison with observations of MCG-6- 
30-15, whereas the new model improves this aspect. 

We also made a comparison of the predictions by the 
new version of the code for the point-to-point variability. 
In this case, the cloud survival decreases the variability 
amplitude since, for the adopted model parameters, about 
half of the clouds survive from one time bin to another and 




Fig. 3. Fractional variability amplitude for a 
Schwarzschild black hole (upper panel) and a maxi- 
mally rotating Kerr black hole (lower panel). Dotted 
lines: uniform, isotropically radiating spots generated 
independently in each time bin. Dashed lines: uniform, 
isotropically radiating spots surviving between con- 
secutive time bins. Continuous lines: present model of 
non-uniform, non-isotropic spots surviving between con- 
secutive time bins. Other model parameters: M = IO^Mq, 
Lx = 10^4 erg s'^, i = 30°, i?out = 50 i?g, Aad = 3, 

7rad - 0, Jrad = 0, iufco = 10^ S, T^bs - 10^ s, iV = 300; 

and nnioan — 30 in the upper panel, while n„ican = 100 
in the lower panel. The inner radius is equal to the 
marginally stable orbit adapted to the given choice of the 
black hole spin. No primary contribution is assumed. 



do not contribute much to the point-to-point variability. 
Some effect is also seen in the shape of Fpp{E) (see Fig.0J. 

Furthermore, we calculated Fyaj{E) and F-pp{E), tak- 
ing into account the contribution of the primary emis- 
sion from the flare, with a relative level, as expected for 
isotropic flare emission. In these cases, the energy depen- 
dence of the variability is weak and comparable to the 
level obtained by Czerny et al. (2004). Although in the 
new approach the energy dependence in the spot emission 
variability is stronger, the spot spectrum contributes less 
to the total (primary plus spot) spectrum due to the less- 
ionized reflecting medium. The two effects roughly com- 
pensate for each other and, for an isotropic primary, the 
variations in _F'vai(£') are again around a mean variability 
level of about 2%. 
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Fig. 4. Point-to-point fractional variability amplitude for 
a Schwarzschild black hole (upper panel) and a maxi- 
mally rotating Kerr black hole (lower panel). Dotted lines: 
uniform, isotropically radiating spots generated indepen- 
dently in each time bin. Dashed lines: uniform, isotropi- 
cally radiating spots surviving between consecutive time 
bins. Continuous lines: present model of non- uniform, non- 
isotropic spots surviving between consecutive time bins. 
The other parameters of the models are as for Fig. O 



3.2. The dependence of the fractional variability on 
the model parameters 

The mean number nmoan of flares mostly affects the overall 
normalization of Fvar, as does the flare life-time distribu- 
tion tiife: larger values of tjife require higher values of n,„oan 
to achieve the same overall variability level. The impor- 
tant ruling parameter for the variability in this context is 
the frequency of flare onsets and endings, which directly 
determines the fluctuations of the lightcurve. If nmoan is 
high, more onsets and endings per unit time happen and 
the averaging effect reduces the variability. The lifetime of 
a single flare acts in the opposite direction. The longer an 
individual flare lives, the less flares set on or end per unit 
time. 

Suppression of i^var at ~ 6.5 keV appears if the radi- 
ation is strongly concentrated toward the center. We see 
such a trend in models with S^e^d — 0, if the spot bright- 
ness radially decreases as the standard disk dissipation 
rate (/3rad = 3) and the flares are distributed uniformly 
across the disk surface (7i.ad = 0). We show an example 
of this behavior in Fig. 13 (continuous line). If /?i-ad is con- 
siderably larger (e.g. /S^ad ~ 4.5), the same pattern is seen 
even for flare life times scaling with the local Keplerian 
value ((5rad — 1.5). This behavior can be qualitatively un- 
derstood from the following considerations: the Ka line 



emitted from the innermost spots is relativistically broad- 
ened and Doppler-shifted according to the local orbital 
velocity vector. Since the life-time of a flare close to the 
marginally stable orbit is, in general, a significant fraction 
of the orbital period, the line emitted by one spot will 
shift smoothly in energy around the line centroid. As the 
line is also strongly relativistically broadened, its shape is 
rather flat around the centroid, which hence exhibits less 
variability than do the much steeper wings. In addition 
to that, there is a contribution to the iron Ka line from 
outer parts of the disk, and its persistent presence further 
supports the drop in i^var around 6.5 keV. The features 
in the variability spectrum are thus mainly determined by 
the relatively few innermost hot spots for which the repro- 
cessed spectrum is more strongly modifled by the general 
relativistic effects. 

The parameters 7rad and (5rad act predominantly in the 
same way, since having more flares at larger distances is 
roughly equivalent to flares lasting longer at large radii. 
However, the effect of the two parameters is not strictly 
identical. The situation is similar with /3i-ad and 7rad: a de- 
crease in /3rad has a similar effect as a decrease in 7i-ad ■ We 
should add here that, to a certain extent, the timescale 
tiiic also influences the energy dependence of Fvar- If the 
timescale is too short, the energy dependence in the vari- 
ability is smeared, which is again due to the large number 
of onsets and endings of the hot spots. 

In summary, the most important property of a given 
parameter set is the sum of the parameters, 2 + 7iad + 
(^rad— /3rad, which represents the index of an effective radial 
distribution of energy generation. If this value is negative, 
the model is likely to show a suppression of the variability 
around 6 keV. If the value is positive, the model is likely 
to show complex behavior around 6 keV in -FVar, as seen in 
Fig. 13 and an enhanced variability at '^ 6.5 keV in i^pp, 
as seen in Fig. [7| However, to some extent the speciflc 
individual values of all parameters are also important. 

As mentioned above, the position of the features in 
Fvar and their shape strongly depend on the relativistic 
effects modifying the spectrum of the inner hot spots. The 
observed features thus constrain both the inclination angle 
i of the disk and the Kerr black hole parameter, a. A larger 
inclination tends to move the features toward higher en- 
ergies, while a faster black hole rotation produces slightly 
broader and much deeper features. This is due to the effect 
that the marginally stable orbit of a rotating black hole is 
smaller, and thus the hot spots can exist at a closer dis- 
tance to the event horizon. Therefore, if any features are 
present in Fvar, they can be used effectively to constrain 
a and i. 

3.3. Comparison to the special case of MCG-6-30-15 
3.3.1. Observational data and fixed model parameters 

Two long X-ray exposures were obtained with XMM- 
Newton for the Seyfert-1 galaxy MCG-6-30-15: one obser- 
vation lasting 95 ks was performed in June 2000 (Wilms 
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Fig. 5. Fractional variability calculated from several 
lightcurves calculated for the same model, with the for- 
mal errors displayed for one of the curves. Upper panel: 

rimean = 125, /3rad = 3; lowCr panel: nmoan = 30, /^rad = 2. 

Other model parameters: AI — IQ'^Mq, a — 0.95, Lx = 



> X 10'*^ erg s 



30°, R. 



50r„ 



0,^r. 



1.5, 



lifco 



2 X 10* s, T 



obs 



out ^^' g: Trad ^i ^rad 

6146 s, and N = 50. No primary 



contribution is assumed. 

et al. 2001), and the second one lasting about 325 ks was 
done in July 2001 (Fabian et al. 2002). The time variabil- 
ity in the data from 2000 was analyzed in detail by Ponti 
et al. (2004), while the time variability in the data from 
2001 was studied by Vaughan & Fabian (2004). 

We attempt to reproduce the shape of the fractional 
variability amplitude obtained by Ponti et al. (2004; see 
their Fig. 2 and Fig. |H1 in this paper). We also derive 
the point-to-point fractional variability amplitude for the 
same data sequence. In our simulations, we adopt a fixed 
value Tobs = 6146 s, as was done for the computations of 
Fvar during the data analysis by Ponti et al. (2004). The 
number of time bins in the i^var calculations is N = 16, 
which is appropriate for the data length. Fpp is calculated 
from the same data using 1000 s bins, so in modeling Fpp, 
we take Tobs = 1000 s, and the number of points increases 
to TV = 100 to represent the data length correctly. We 
adopt lO^A/0 for the mass and 8 x 10^^ erg s~^ for the 
X-ray luminosity of the source. We fix the outer disk ra- 
dius somewhat arbitrarily at i?out = 50 Rg. In most cases, 
we assume i = 30°, a Kerr parameter of a = 0.95, and 

7rad = 0. 

3.3.2. Intrinsic randomness of the model 

The most noticeable property of the model is the consid- 
erable randomness in the results for some sets of model 



Fig. 6. Fractional variability calculated independently for 
three separate parts of the 300 ks lightcurve of MCG-6-30- 
15 (XMM-Newton data, 2001, thin points). The set with 
filled circles represents the result obtained from the 100 
ks hghtcurve (XMM-Newton data, 2000). The length of a 
time bin equals 6146 s. 

parameters and the relative stability for others. The ef- 
fect is shown in Fig. |5| We consider two exemplary sets 
of parameters. In the first set (upper panel in Fig. EJ, the 
spot brightness decreased with the radius as /?rad = 3, in 
the second set with /3rad = 2. Again, we choose the mean 
number of clouds in each model in such a way as to have 
roughly the same level of variability. For a specific choice 
of the other parameters (the same for the two models; see 
caption of Fig. El, we have to assume nmoan = 125 for the 
first model and nmcan = 30 for the second one. 

For both parameter sets, we obtained five 300-ks-long 
time sequences and calculated Fvar for each of the se- 
quences independently. For clarity, the expected statistical 
one sigma errors are marked for one of the curves only. The 
curves in the upper panel are roughly similar, and they do 
not differ more from each other than expected from sim- 
ple error analysis. However, the curves in the lower panel 
differ rather widely. This demonstrates that, for such a 
set of parameters, a single lightcurve of a standard length 
and with a standard error description cannot be used to 
estimate whether the model well represents any specific 
data. The effect is still somewhat stronger for our data 
covering only 100 ks instead of 300 ks, as was adopted in 
these simulations. 

The discrepancy is due to the different total mean 
numbers of flares required by the two data sets. Since 
the flare duration distribution was the same in the two 
models the total (time-integrated) number of flares in the 
lightcurve is 726 for the first set of parameters and only 
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Fig. 7. Point-to-point fractional variability calculated for 
the same models and displayed in the same manner as the 
fractional variability shown in Fig. [S] 

174 for the second set. The standard determination of the 
statistical error (see Eqn. A. 2 of Ponti et al. 2004 and 
Eqn. in this paper) of course does not include this ef- 
fect. This result supports the conclusion of Vaughan et al. 
(2003), based on Monte Carlo simulations of lightcurves 
with an assumed power spectrum, that the overall nor- 
malization of the variance is rather unreliable if based on 
short observations. However, our energy-dependent calcu- 
lations indicate that although the overall level of the i^var 
changes strongly for some sets of parameters, the shape 
of -Fvar varies less from sequence to sequence. The energy 
pattern is roughly preserved and can be used to estimate 
the source properties. 

Since our conclusion is based on a specific flare/spot 
model, we compare the trend to the sampling effects in real 
data: in Fig. IHl we show three F^ar (E) plots obtained by 
the division of a single 300 ks long lightcurve of MCG- 
6-30-15 (XMM-Newton observation performed in 2000, 
see Fabian et al. 2002) into three equal parts. Fvar(£') 
was calculated for each of the three parts independently. 
The overall change in the variability level between the 
three subsequences is in good agreement with the expected 
power red leak for the source. Table 1 of Vaughan et al. 
(2003) for a source with a PSD slope of 1.0 below the break 
and 2.0 above the break, for 20 data points, predicts the 
scatter between -0.71 and 0.45 in logarithm of the vari- 
ance, so between -0.36 to 0.23 in logarithm of dispersion. 
This means that the average value of 0.2 is predicted to 
fluctuate between 0.09 and 0.33, as observed (see Fig. EI). 
However, most of the variability is in the overall normal- 
ization of _F'var(£'); the shape actually varies much less, al- 
though some changes in energy dependence are also seen. 



Fig. 8. Point-to-point fractional variability calculated in- 
dependently for three separate parts of the 100 ks 
lightcurve of MCG-6-30-15 (XMM-Newton data, 2000, 
thin points). The set with filled circles represents the mean 
result obtained directly from the whole lightcurve. The 
length of a time bin equals 1000 s. 

The variability is always suppressed just above 6 keV, but 
one of the data sets does not show the complex structure 
in 4 - 6 keV band. 

The structure of the F^^^{E) spectrum of MCG-6-30- 
15 shown in Fig. is quite typical for Seyfert galaxies. 
The suppression around 6 keV is often found to be even 
more pronounced than for this particular object (see, e.g., 
Markowitz, Edelson, & Vaughan, 2003, for some exam- 
ples). Another important feature of the rms spectra is 
an increase of the variability in the soft X-ray band. An 
interesting interpretation of this feature is proposed by 
Chevallier et al. (2006). According to this idea, the en- 
hanced soft X-ray variability is due to the response of 
the warm absorber to the variations of the incident flux. 
Absorption is most efficient around ~ 1 kcV and the 
change of the ionization state of the warm absorber in 
the response to the variable flux amplifies the overall vari- 
ability level in this spectral band. This leads to an increase 
in the simulated fractional variability amplitude at about 
that energy. 

We also analyze the scatter in point-to-point fractional 
variability amplitude predicted by the model. The result 
is shown in Fig. [7| It is interesting to note that the results 
are complementary to the results for the simple fractional 
variability. While in Fig. [3 there is much more scatter in 
the upper panel, in Fig.[7|the trend is reversed. This com- 
plementarity is connected with the presence of relatively 
more (or less) power in a given model at longer timescales. 
Also, in this case, the model predicts much less variabil- 
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ity in the shape than in the overall normalization. In the 
observed Fpp, calculated separately for three parts of the 
lightcurve, we notice considerable variations in the nor- 
malization, although lower than for i^var (see Fig. IS)). We 
see some variations of the shape as well, but the basic 
pattern remains the same. There is a decrease of the vari- 
ability at ~ 6 keV, which is present in all three sequences, 
although in one of them it is shifted from 6.6 keV to 6.0 
keV. It is difficult to discuss the significance of such vari- 
ations since the errors are large, although the number of 
bins was relatively high due to the choice of short time 
bin (1000 s). 

This analysis shows that the observational errors are 
actually much larger than indicated by a simple statistical 
analysis of a given single lightcurve, as is expected from 
the power leaking (Vaughan et al. 2003). However, the 
errors in the shape are not as large as the overall change of 
fractional variability amplitude from one to another time 
sequence, due to the long timescales present in the system. 
Of course, an increase of the duration of the observation 
time would give less scatter. But then, such requirements 
may be difficult to meet. The model from the lower panel 
of Fig. gives results for i^var that are almost within the 
statistical errors of a single lightcurve, if the number of 
time bins is increased to iV = 300. This corresponds to an 
observation of more than 20 days. 



3.3.3. Best representations of Fpp(£') 

In Fig.|51 we show the point-to-point fractional variability 
amplitude of MCG-6-30-15 and the best representation of 
the short timescale variability from our model. The model 
well represents the suppression of the variability near the 
core of the iron line, in the 3-10 keV energy band. The 
position of the minimum in the theoretical plot is basically 
determined by the inclination angle, and apparently, the 
value of 30° adopted in most simulations is appropriate. 
Inclinations smaller than ^ 10° or larger than ^^ 40° are 
excluded. 

The width of the dip is mostly determined by the Kerr 
metric parameter a. Fast rotation of the black hole results 
in a broad feature extending far to the red side of the line 
core, while a slowly rotating black hole gives a narrow fea- 
ture. The value a — 0.95, relatively high, but still signifi- 
cantly lower than the maximally rotating black hole, best 
represents the data. A slower rotation than a ^ 0.8 is not 
consistent with the data. The depth of the feature mostly 
depends on the ratio of the spot spectrum to the primary 
emission going towards an observer. Our best representa- 
tion was obtained assuming no primary contribution, but 
the contribution of the primary at some level is allowed, 
as analyzed in detail in Sect. 13.331 The overall level of the 
variability depends mostly on the mean number of flares, 
and this can always be adjusted if the energy-dependent 
shape of the i^var is appropriate. However, the same model 
parameters do not represent the Fvar well. 



0.4 




Fig. 9. The best model of Fpp (continuous line) with the 
corresponding data (lower set of points) . The upper data 
set represents i^var and the modeling result for the same 
parameter set: /^rad = 4, 7rad = 0, ^rad = 1-5, tiifc,, = 
2 X 10^ s, and nmcan = 750. The other parameters are at 
their standard values (see caption of Fig. ^ . The model 
is shown unbinned for better clarity. 



3.3.4. Best representations of F^^riE) 

In Fig.^1 we show the best representation of the standard 
Fvar{E) for the source. It roughly reproduces the observed 
'wiggles' of the i^var obtained from the whole lightcurve. 
However, this time the same parameters do not reproduce 
Fpp correctly. Since there is a considerable randomness 
in the modeling process, we repeated the computations 
several times for a different random seed of the Monte- 
Carlo routines. But in all such realizations, the Fpp shows 
an excess at 6.4 keV instead of a dip. A discrepancy with 
the data therefore remains. 

To decide which set of parameters better represents the 
source properties, we plot the spectra and the lightcurves 
from the two models (see Fig. Ill|) . The spectra of both 
models are quite comparable and do not help to differen- 
tiate between the candidate models. The lightcurves, on 
the other hand, have a rather different character visually 
(see Fig. ^J . The upper curve shows less short time-scale 
flickering and seems to be quasi-periodic, while the lower 
curve looks more similar to the source X-ray lightcurve 
(see Fig. 1 in Ponti et al. 2004). 

Therefore, we consider the set of parameters used for 
Fig. |51 as a more satisfactory description of the variability 
in MCG-6-30-15. 
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Fig. 10. The best model of Fvar (continuous line) with the 
corresponding data (upper set of points) . The lower data 
set represents Fpp and the modeling result for the same 
parameters: /3rad = 3, 7rad = 0, 5rad = 1-5, iufeo = 2 x 10*^ 
s, and rimean = 100. The other parameters are at their 
standard values (see caption of Fig.|31). The model is shown 
unbinned for better clarity. 



3.3.5. Contribution of the primary emission 

To directly compute observable characteristics, we need 
to combine the primary continuum with the reflected 
component. The resulting signal depends on the mutual 
proportion of the two components, and also on the en- 
ergy range and time resolution of the observation. This 
means that the degree of anisotropy of the primary flare 
should be taken into account. Our modeling shows that, if 
the primary power law emission strongly dominates over 
the reprocessed component, and if the slope of the pri- 
mary emission does not change with time, there is no 
energy-dependence in the fractional variability amplitude. 
Moderate contribution of the primary results in flattening 
of the energy dependence without changing the overall 
variability amplitude. We note, though, that this is only 
true for the case of a well-defined correlation between the 
primary and the reprocessed component. This correlation 
implies that the flare geometry shown in Fig.^is the same 
for all flares distributed over the disk. In reality, these 
things are likely to be more complicated. If the primary 
emission and the reprocessing are not perfectly correlated, 
the overall level of the variability goes down when the pri- 
mary component becomes more significant. 

We calculate examples of the model shown in Fig. |51 
with different levels of primary contribution and assuming 
that both components are well-correlated. We also calcu- 
late the equivalent width of the iron line in the resulting 



Fig. 11. The spectra at various moments computed from 
the models shown in Fig. IHl (lower panel) and Fie. 1101 fup- 
per panel) , with the mean spectra over-plotted by a thick 
line. 

Table 1. Equivalent widths (EW) of the iron line complex 
and the ratio -F'pp(6 keV)/i^pp(6.6 keV) as a function of the 
amount of reflection. The observed i^pp(6 keV)/Fpp(6.6 
keV) ratio is 1.19, and the model results were binned into 
the same energy channels as the data. 
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spectra. These results are shown in Table Q] The model 
without any primary contribution (anisotropy parameter 
TZ = oo) best represents the decrease of the variability at 
~ 6.6 keV, but the EW of the iron line complex in the 
predicted mean spectra is very large, as can be seen in 
Fig. ^] The reflection amplitude 7?, ~ 4 is still roughly 
acceptable from the point of view of the variability. The 
actual measurements of the line EW in MCG-6-30-15 gave 
values ranging from 800 eV to 250 eV, depending on the 
assumption about the underlying continuum (Wilms et 
al. 2001), and Reynolds et al. (2004) obtained the value of 
500^200 "^^i ^'^ ^he line intensity is not strongly overesti- 
mated. The model with TZ ~ 2 seems to have too low of an 
energy dependence in the fractional variability amplitude, 
and for 7?, = 1 it underestimates the line intensity as well. 
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Fig. 12. The lightcurves of the models shown in Fig. 
(lower panel) and Fig. 1101 (upper panel) for narrow energy 
bands around 4.4 keV, 5.2 keV, 6 keV, 6.4 keV, 6.8 keV, 
and 8.8 keV, and with the resolution used for the Fpp 
computations. 



The amount of reflection can be estimated from the 
strength of the Compton hump seen in hard X-ray data. 
The analysis of Beppo-SAX data by Ballantyne et al. 
(2003) showed that the reflection factor determined from 
PEXRAV was R = 2.7^09 ^'^^ solar abundances. Assuming 
higher iron abundance than solar improved the fit and only 
gave a lower limit to the reflection factor (i? > 2.6) so the 
spectrum can actually be reflection-dominated, as advo- 
cated for NLSl (including MCG-6-30-15) by Fabian et al. 
(2002). A purely reflection-dominated spectrum {R > 3) 
for MCG-6-30-15 was also considered by Taylor, Uttley, & 
McHardy (2003). Therefore, from the point of view of the 
primary contribution, the favored model is satisfactory. 

4. Discussion 

The X-ray spectra of a number of AGN contain a reflection 
feature which is best modeled as a relativistically smeared 
iron line (for a review, see, e.g., Reynolds & Nowak 2003). 
If the source of the irradiation consists of multiple local- 
ized flares, it is natural to expect the variability of a single 
spot luminosity (e.g. Abramowicz et al. 1991; Bao 1992; 
Karas et al. 1992), as well as of the line shape (Ruszkowski 
2000; Nayakshin & Kazanas 2001; Turner et al. 2002; 
Iwasawa, Miniutti, & Fabian 2004). 

In the present paper, we studied in detail the time- 
dependent spectra from a multi-flare model, paying par- 
ticular attention to the energy-dependent fractional vari- 
ability amplitude as a diagnostic tool. We used both 
the standard fractional variability amplitude, Fvar, and 



the point-to-point fractional variability amplitude, Fpp. 
We showed that the model explains a decrease in the 
variability level at ~ 6.5 keV derived from the rms- 
spectra of Seyfert galaxies, such as MCG-6-30-15 (Inoue & 
Matsumoto, 2001; Vaughan & Fabian, 2004; Ponti et al., 
2004) . Similar behavior was observed in Mrk 766 (Pounds 
et al., 2003), 3C 120, Akn 564, IC 4329a, 3C 390.3, 
FairaU 9, NGC 3516, NGC 3783, NGC 4151, NGC 5548, 
or NGC 7469 (Markowitz, Edelson, & Vaughan, 2003). 

The overall variability level, predicted by the multi- 
flare model and expressed in terms of i^var and Fpp, de- 
pends strongly on the adopted number of flare-spots. By 
adjusting this number, the model is capable of reproduc- 
ing the overall normalization of observed rms-spectra as 
we demonstrated for the case of MCG-6-30-15. However, 
the form of the energy dependence of Fvar and Fpp is more 
difficult to adjust and provides a strong diagnostic tool ei- 
ther for proof or for falsification of our model. 

The prediction of the multi-flare model for the shape 
of Fvar and Fpp depends predominantly on the radial pro- 
file of the energy dissipation characterizing a given set of 
parameters (see Sect. 13. 2|) . If the energy dissipation in the 
form of flares weighted within the emission area decreases 
with the radius, the model shows a suppression of the vari- 
ability at 6.6 keV. 

The present data predict both Fvar and Fpp with sig- 
nificant error, mostly due to the presence of the timescales 
much longer than a typical lightcurve duration (e.g. 
Vaughan et al. 2003). The errors in the shape of the energy 
dependence seem to be lower, but still rather difficult to 
determine observationally. Therefore, we did not provide 
formal fits to the observational data at the present stage. 
Instead, we favored one of the studied sets of parameters 
at the basis of semi-quantitative analysis. 



4.1. Best model properties for MCG-6-30-15 

The parameters of the model which we considered satis- 
factory (see Fig.|21l are interesting. The distribution of the 
timescales was Keplerian, and the requested normalization 
leads to a timescale range from 1910 s in the innermost 
part to 8.7 x lO'' s in the outermost part. The spot size 
obtained was i?x = 1-96 Rg, which is not too large in the 
context of the underlying model assumptions. It corre- 
sponds to a height of the flare source above the disk equal 
to 1.13 Rg. The mean number of coexisting flares is quite 
large (750), but the variability is nevertheless observed 
to be that large since most of the luminosity comes only 
from the few innermost flares. This is due to the rather 
steep increase of the flare luminosity toward the disk cen- 
ter (/3rad = 4). Interestingly, the obtained mean number 
of flares is not unlike the numbers of clouds obtained in 
the obscurational variability model of Abrassart & Czerny 
(2000). 

The steep flare emissivity mentioned above is indeed 
higher than might be expected just from the shape of the 
gravitational potential. However, modeling of the mean 
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shape of the broad iron Une frequently required an en- 
hanced emissivity profile (e.g. Wilms et al. 2001). It was 
also predicted by a theoretical corona model of Kawanaka 
et al. (2005). 

Although the mean number of the coexisting flares is 
large, we see occasional large 'single' flares in the theoret- 
ical lightcurve. Such resolved single flares are already seen 
occasionally in the present data (e.g. Ponti et al. 2004 for 
MCG-6-30-15). Turner et al. (2005) identified a single or- 
biting fiare in Mrk 766 on the basis of the spectral analysis, 
and Porquct et al. (2004) also considered a localized flare 
as a plausible explanation of the ESO 113-GOlO spectrum. 

4.2. Model approximations 

There are still certain inconsistencies present in our model, 
mainly in the computation of the local spectra. The calcu- 
lations were made for a Schwarzschild black hole of 10^ so- 
lar masses. However, in computing the relativistic effects, 
we also use them for other values of M . For MCG-6-30-15, 
a smaller mass is more appropriate. In our computations 
of the relativistic effects, we assumed l(f Mq for the mass 
of the central black hole, and McHardy et al. (2005) sug- 
gested an even smaller value ^ (3 — 6) x lO^M©). However, 
we do not expect any strong direct dependence of the local 
spectrum on the black hole mass. 

The local computations were performed at a single ra- 
dius (18i?g), but we checked that the local spectra cal- 
culated at 7i?g do not differ strongly (Goosmann et al. 
2006) when computed for the same flare-to-disk-flux ra- 
tio below the flare. Therefore, for models calculated with 
this ratio preserved (i.e., models with /Jrad = 3) we do not 
think we introduce much of an error. However, for models 
with different /3rad the local spectrum may be different, so 
a range of local models should actually be calculated to 
incorporate this effect. 

The assumptions about the primary radiation in our 
model are rather simple. The distance between the flare 
source and the hot spot is assumed to be negligibly small 
so that the relativistic effects acting on the incident radia- 
tion can be neglected. In Goosmann et al. (2006), we show 
that this is a good approximation for Schwarzschild black 
holes. For Kerr black holes it works fine if the flares do 
not get closer to the black hole than ^ 3i?g. The primary 
source is assumed to turn on and off instantaneously and, 
consequently, so arc the hot spots. We thereby assume that 
the time for the actual reprocessing of the radiation in the 
disk medium is negligible. The local lightcurve of the spot 
emission hence has a simple box shape. We keep this as- 
sumption, although it is probably too simple, since not 
much is known yet about the actual production of the pri- 
mary radiation and its time evolution during a magnetic 
reconnection event. In later developments of the model, 
more detailed assumptions about the nature of the flare 
source should be included. Also the ratio of the primary 
flux to the disk flux is somehow arbitrarily set to 144 at 
the spot center. In future work, this flux ratio should be 



treated as a free parameter, but of course this requires a 
large grid of calculations for the local spectrum. 

4.3. Power spectrum 

The power spectrum derived by McHardy et al. (2005) for 
MCG-6-30-15 shows a bending, characterized by the break 
frequency 3.2 x 10~^ Hz, and the asymptotic slopes of -0.8 
and -1.74 if the broken power law description is used. It is 
not clear that our model is able to reproduce its properties. 
We attempted to calculate the power spectrum from our 
theoretical lightcurves, but the result was consistent with 
white noise. This may be partially due to the fact that: (i) 
the lightcurves contained only a limited number of points 
(100), and (ii) the studied dynamical range of the model 
(-Rout = 50 i?g) is not yet fully satisfactory. 

On the other hand, our model in its present version 
may not contain enough power in the long frequency 
range. Generally, more complex models involving corre- 
lated avalanches of flares are used to model the power 
spectra in AGN and galactic sources (e.g. Poutanen & 
Fabian 1999; Zycki 2003; Zycki & Niedzwiecki 2005). A 
successful model accounting for the power spectrum has 
to explain why there is a lot of power in the low frequency 
part of the power spectrum, while most of the energy dis- 
sipation takes place in the innermost part where the ex- 
pected timescales are short. 

Two essentially similar models were proposed which 
can explain the presence of the long timescales success- 
fully: 'self-organized criticality' (Mineshige et al. 1994), 
and 'traveling perturbations', (Lyubarskii 1997). Both 
conclude that the disk accretion rate is locally perturbed, 
and this perturbation propagates inward, so the accre- 
tion rate at any disk radius contains a history of long 
timescale perturbations created at larger radii. In partic- 
ular, the model of Lyubarskii (1997) looks natural within 
the frame of the magneto- rotational instability mechanism 
of viscosity, and it is particularly appropriate for a power 
spectrum with a single break like that of MCG-6-30-15. 

Our flare model, at present, does not include this kind 
of perturbation memory. However, in principle, this kind 
of non-local perturbation of the disk accretion rate can be 
incorporated into the model. If it is combined with some 
assumptions of scaling between the local momentary ac- 
cretion rate and a flare luminosity, it will probably provide 
the missing power at long timescales. 

4.4. Ottier models of X-ray variability 

We have found that our multi-flare model with a time- 
independent shape of the energy spectrum of the flare 
emission provides an interesting interpretation of the drop 
in the variability around 6.5 keV. However, there are other 
models under discussion. 

The energy-dependent form of the variability was also 
modeled as a variable primary emission with constant 
reflection (e.g. Taylor, Uttley, & McHardy, 2003; Shih, 
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Iwasawa, & Fabian, 2002), or as a result of the pivoting of 
the primary power law emission (e.g. Markowitz, Edelson, 
& Vaughan, 2003). 

However, in this first class of models it is difficult to 
explain why the reflection component should be constant 
since either it forms close to the black hole (and the source 
of the primary emission) and should respond to changes 
of the primary, or it forms far away, then it is difficult 
to explain why it is relativistically broadened. An attrac- 
tive way out of this paradox was proposed by Martocchia, 
Karas, & Matt (2000); Miniutti & Fabian (2004); and 
Miniutti ct al. (2005), who describe a model where the 
change in the primary intensity is predominantly due to 
the minor changes in the position of the source of the pri- 
mary radiation at the disk rotation axis, which gave effects 
strongly enhanced due to the general relativity. Pivoting 
is also quite natural, if the flare evolution (cooling, expan- 
sion) is taken into account. However, in this case the drop 
in the variability level can be at any position, and not just 
around the iron line complex. 

There are also entirely different models of the X-ray 
variability of AGN, for example the presence of the den- 
sity inhomogeneities in the disk atmosphere resulting from 
the photon bubble instability (Ballantyne et al. 2005), de- 
velopments of the temporary spiral structures in the ac- 
cretion disk (Fukumura & Tsuruta 2005), variations in 
the properties of the warm absorber (e.g. Sako et al. 2002; 
Done & Gicrlinski, 2006, private communication), or vari- 
ations in the the inner disk torque (Garofalo & Reynolds 
2005). Only more detailed studies of model predictions 
and comparison with the data in more than just one as- 
pect may bring a reliable answer to the true nature of the 
X-ray variability. 
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